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Abstract 

Stark broadened lines of the hydrogen Brackett series are computed for the conditions 

of stellar atmospheres and circumstellar envelopes. The computation is performed within 

l^j the Model Microfield Method, which includes the ion dynamic effects and makes the bridge 

between the impact limit at low density and the static limit at high density and in the line 

> wings. The computation gives the area normalized line shape, from the line core up to the 

static line wings. 
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1. Introduction 



Hydrogen is the most abundant element in the Universe. Its broad lines give noticeable 
features in the spectra of stellar atmospheres [U [21 El H]- These lines are very sensitive to 
the interaction between hydrogen radiating atoms and the surrounding charges, electrons 
and ions (mostly protons), which is connected to the random electric field generated by 
these charges. The electric field has two components, with different time scales: the rapidly 
varying electronic field and the slowly varying ionic electric field. The net field induces a 
strong mixing of the atomic states with the same principal quantum number, from which 
the Stark broadening originates. Astrophysical applications need to know the full line shape, 
from line centre up to the line wings, for a wide range of plasma conditions (i.e. electron 
density N e and temperature T). 

Whereas hydrogen line shapes have been widely used in stellar physics for the determi- 
nation of gravity and/or temperature in the visible and UV ranges, the recent instrumental 
developments on various telescopes, such as AMBER [5| and CRIRES [6] on VLTI and VLT, 
require the availability of precise hydrogen lines in the infrared (i.e. between 1 and 5 //m). 
The latter are badly known [?]. 

This paper thus aims at providing a coherent description of the line shapes of Brackett 
a, P, 7, which connects the levels of principal quantum numbers n equal to 4, for the lower 
state of the transition, and n! equal to 5, 6 and 7, for the upper state, and have central 
wavelengths of 4.05, 2.63, and 2.12 /im. The choosen plasma conditions are relevant to 
stellar photospheres and circumstellar environments (electron densities between 10 10 and 
10 19 cm~ 3 and temperatures between 10 3 and 10 7 K). The plasma charges are assumed to be 
electrons and protons. This is a standard approximation for this type of study, although 
some improvements may be possible by including the effect of ionization of He (which is, by 
number of atoms, only 10% as abundant as hydrogen) and traces of heavier elements. 

We suppose that plasma collective effects are included in the Coulomb interactions be- 
tween the hydrogen bound electron and the plasma charges by Debye screening. This requires 
that the Debye length \ d , is larger than the mean distance r between the protons, or that 
the parameter a = r /\ d is smaller than unity. Our tables will thus be limited to the values 
of temperatures and densities satisfying the condition: 



which is fulfilled for standard stellar atmospheres. 

We assume also that the proton/electron density is small enough to ensure that each line 
4-n' (n'= 5, 6 or 7) remains distinguishable from the subsequent line 4-(n'+l) of the Brackett 
series. Using the Inglis-Teller |8] criterium, the upper limit to the electron density N e is thus 



This respectively gives log(A^ ei7naa; ) = 16.75,16.16 and 15.66 for the n' (= 5, 6 and 7) 
values considered in this paper. However, in order to allow interpolation within the tables 
for astrophysical applications, the line shapes for higher values of N e , up to a decade, have 
been computed. 



N e < 2.6 10 9 T 3 (cm" 3 , K), 



(1) 



log(N e , max ) = 22- 7.51og(n') (cm" 3 ). 



(2) 
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Different methods can be used to generate high quality spectral line shapes for hydrogen 
lines perturbed by protons and electrons: Molecular Dynamics for describing precisely the 
ion dynamics effects in the line center |9l[T0], quantum theory for the electron contribution to 
the line wings [11], short range H-H+ molecular interactions also for the line wings, leading 
to the apparition of quasi-molecular satellites, which are observed in the atmospheres of 
white dwarfs [12]. However, they are limited either to part of the profile, or to restricted 
plasmas conditions, or to simple lines like Lyman or Balmer lines. Thus, as for astrophysical 
purposes, the tabulations go from the line centre up to the line wings, it is necessary to find 
a compromise between accuracy and description of the whole profile. 

In this context, the tabulations of Vidal et al.[T3], using Unified Theory for the electrons 
and static approximation for the ions, have been used for stellar atmospheres, despite the 
intrinsic lack of accuracy in the line center due to the neglect of ion dynamics effects. The 
tabulations of Stehle et al., for the Lyman , Balmer [TU [15] and Paschen lines [E], using 
Model Microfield Method, which was initially developed by Brissaud et Frisch [TTlttE] brought 
an important improvement by taking into account the ion dynamics effects. They are now 
used for atmospheres and for the computation of radiative diffusion processes in the radiative 
stellar envelopes [H]. In the case of partly ionized plasmas, like for the atmospheres of cool 
stars, the contribution of neutral broadening by hydrogen has to be included in the line 
shape, especially in the line wings [201121]. This effect will be neglected in the following. 

In this paper we shall present the Stark broadened profiles of Brackett lines. They will 
be computed within the formalism of Model Microfield Method, hereafter denoted by MMM. 
We shall neglect the fine structure effects, which play a role, at low density, in the core of 
the lines with low n quantum mumbers, like Lya or Ha [22]). 

For Bra , the profiles are computed at densities logio (N e (cm -3 )) ranging from 10 to 18.5, 
by step of 0.5. For each density, the profiles are computed at temperatures equal to 1000, 
2500, 5000, 10000, 19550, 39810, 79810, 158500, 316200, 63100 and 1259000 K, assuming 
that condition [j] is satisfied. For Br/3 , the computation stops at logio (N e (cm -3 )) = 18 and 
for Br7 at 17.5, as explained previously. 

2. Method 

The broadening of spectral lines results from the interactions between the radiating hydrogen 
atom and the free ions and protons. These two contributions can be described in terms of 
interaction potentials, with the corresponding electronic and ionic plasma microfields F e / and 
F ion . Neglecting quadrupolar and other contributions that play a role at high densities [23j, 
the dipolar potential of interaction between the bound electron and the microfields may be 
written as: 

V(t) = -d.(F el (t) + F ton (t)). (3) 

The spectral line profile I(uS) (with area normalized to unity) is thus defined in the Liou- 
ville space [2UI2S], spanned by the states \i, f 3> (which stands for |nj, Z i; nf, If, rrif 
as 
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J ^ = ^Td d^^ 6 ^ <, '•/• <1 ' : • /• < U H >el,ion;if;i'f>, (4) 

f-f lJ ' l 'f i,f,i',f 

where < U(uj) > e i,i on is the Fourier transform of the evolution operator of radiating 
Hydrogen atom in the Liouville space, averaged over the realizations of the stochastic dynamic 
electronic and ionic microfields F e i and F ion . The term djj.dj/j/ is the product of dipole 
operator elements between initial low states (denoted by i, i') and final upper states (denoted 
by /, /'), of the hydrogen bound electron. As fine structure and inelastic effects are neglected, 
one has Ef — E{ = E' f — E[ = hw . 

The two microfields are stochastic processes. It is thus possible to define two distribution 
functions -P(F) [26, 27J, respectively associated to the slowly varying ionic and rapidly varying 
electronic microfields. In order to take into account the dynamic effects of these microfields, 
a model for the dynamical statistics of field fluctuations is necessary. In MMM, the microfield 
fluctuations are handled with a statistical process model, where the microfield (electronic or 
ionic) is assumed to be constant during a given time interval. The microfield then jumps 
instantaneously to another constant value for the next time interval. The jumping times 
are assumed to follow a Poisson law, with a field dependent frequency v(F). The jumping 
frequency v(F), is chosen to reproduce the true field autocorrelation function [17], [JS], [28j. 
This method has been tested against asymptotic impact and quasi-static limits and has been 
proved to lead to very good results for hydrogen [29j and hydrogenic ion lines [28] . 

The method has been already described in Stehle and Hutcheon [15], and we refer the 
reader to this paper for the details. An important point is that it is possible to disentangle 
the contributions of ions and fast electrons by introducing a frequency dependent electronic 
relaxation operator 7 e ;(co>), which is independent from ionic fields, and thus may be computed 
separately. The Fourier transform of the evolution operator, < U(u>) > e i, averaged over all 
the realisations of the electronic fields, may be written as 

< U(w) > el =t[AuoI + t lel (uj)}- 1 , (5) 

where Auj = to — u> Q is the detuning from line center and I is the identity operator in the 
Liouville space. 

Thus, this electron damping is first computed to account for average effect of the electronic 
fields. Then, the static Fourier transform of the evolution operator, < U(oo) > e i,ion, averaged 
over the realisations of the electronic and ionic fields, may be written as : 

< U(cu) > e i, ion = i J P{F lon )dF lon [Aul + d.F lon + ^{u)]' 1 . (6) 

The MMM expression is more complex ( [TT1, [T8] ) than this one, which corresponds to the 
usual Unified Theory with static ionic fields. 

In order to reduce the dimensions of the Liouville space, (16 x 25 = 400 states for Bra, 
for instance), we use the formalism of the reduced Liouville space, which takes advantage of 
the invariance of the different operators, like d.d, by angular average over all the orientations 
of the electric fields and of the fact that the dipole tensor d in equation [3] is of rank 1 (see for 
instance [281 [29]). However, the number of reduced states \rii,li;nf,lf 3> (with \l { —lf \ < 1) 
remains important (i.e. 10 reduced Liouville states of rank 1 for Bra). Thus, we use another 
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approximation, already called « isotropic » approximation in [To\, which uses the diagonal 
form 7gf°(u;), with equal diagonal matrix elements, instead of j e i(uj). This scalar tensor is 
deduced from the pure electronic profile I{oj) e i by the relation, 

I(u) el = --Im[Auj + l yi°(uj)}- 1 , (7) 

7T 

or 

E d i,/- d *',/'7e«;i/;i'/'( W ) 

7$» = ^ VH d . • (8) 

i,f 

The Stark profile is thus obtained after computing the average over electronic field values, 
which gives 7**° (u;), and then the average over the ionic fields, both using the formalism of 
Model Microfield Method. The third step is then the Doppler convolution, which gives the 
final line shape. 



3. Results 

The three Brackett a, j3, 7 line shapes have been studied in a wide range of stellar conditions, 
but we report hereafter only a selection of results, relative to the line center and the line 
wings. 

3.1 Line center 

We present in figures [TJ [2] and [3] the variations of the half- width (HWHM) of Bra, Br/? and 
Br7 lines for different values of the electron density. This quantity is the detuning from the 
line center, at the point where the profile reaches half the maximum value of the line profile 
(which is not necessarily at the line center, as will be discussed below). The figures show the 
half width of the MMM profile with and without Doppler effect and also the value of the 
impact half-width in its own validity range. 

Hence, impact limit has been proved to be valid, for both electron and ion contributions, 
in the line center and at low densities. Moreover, it has been proved that the value of the 
impact width is analytical and that the corresponding profile is Lorentzian in the line center 
|30] . The validity condition of the impact limit is that the half-width value is smaller than 
the ion plasma frequency. Let us take the example of of Bra. At 10 12 cm~ 3 and 2500 K the 
ion plasma frequency is equal to 2.6 10 9 rd.s" 1 , whereas the ion impact contribution to the 
HWHM is equal to 1.8 10 9 rd.s -1 . At 10 14 cm -3 , they are respectively equal to 2.6 10 10 and 
3.9 10 10 . The impact limit should thus be reached gradually as the density decreases below 
10 13 cm" 3 . 

The figure [T] relative to Bra, shows indeed that the half-width of the Stark profile (black), 
tends to converge towards the impact analytical limit (red) at these low densities. However, 
the convolution with the Doppler profile increases the half-width value. As a consequence, 
the half width is dominated by the Doppler broadening at low densities. Similar behaviour 
occurs for the other lines, as may be on figures [2] and [3] for Br/3 and Br7. 
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Figure 1: HWHM of Bra at 2500 K, in A, versus electron density N e in cm 3 (black and 
squares: MMM Stark only; red and circles: analytical impact formula; dashed green: MMM 
Stark profile with Doppler convolution) . 



At higher densities, the ion dynamic effects become smaller in the line center (they are 
negligible in the line wings, as will be seen below), and the line shape departs from the 
Lorentzian shape. We found that, depending on the temperature conditions, the Brackett (3 
line may present a small dip in the line center at moderate densities of 10 15 cm -3 and for the 
largest temperatures as seen in figure |4j This effect was well known in earlier tabulations for 
the Lyman/3 and H/3 lines. The dependence versus the temperature is a consequence of the 
the electron impact broadening, which varies in T -1 / 2 , and which may fill (at low T values) 
or not (at large T values) the central dip. 

3.1 Line shapes 

Another typical behaviour of hydrogen lines is the convergence towards the (Holtsmark) static 
limit in the line wings, which scales as |Aa;| _52 when the line shape is expressed in angular 
frequency units (which is the appropriate unit for the line shape computations). However, 
traditionally, the line intensity was expressed in units of Aa = AX/F , where F (esu) = 
1.25 10 9 (iV e ) 2//3 is the normal Holtsmark field. In these units, the Holtsmark limit is given 
by : 
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Figure 2: 5ame as figure\I\for Br/3. 



J(Aa) 



A /F + Aa 



1/2 



for |Aa| << A /F 



(9) 



|Aa| 5 /2 ^ Ao/F J " |Aa| 5 / 2 
where Aq is the central wavelength, and K a a constant, which depends on the transition 



K a = 1.512 for Bra, 2.401 for Br/? and 2.926 for Br 7 (10) 

This variation in Aa of eq. [9] introduces, at large detunings, a "trivial" asymmetry 
between /(Aa) and /(-Aa) (which does not exist between /(Aa;) and /(— Aa;)). The figures 
5]{7]show the profiles /(Aa) of Bra, Br/3 and Bvy lines at 10 12 cm -3 and various temperatures. 
The profiles, including Doppler effects are reported in red color, and the pure Stark profiles 
in blue. 

The profiles are area normalized, i.e. : 



+oo 



/(Aa) x d(Aa) 



(11) 



As a consequence, broad profiles have small values of /(Aa = 0). For the density consid- 
ered in the figures [5]{7j the Doppler profile dominates the Stark profile in the line center and 
becomes indistinguishable, as expected, in the line wings, where they follow the asymptotic 
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Figure 3: Same as figure\I\for Bry. 

limit of eq. [9j At higher densities (not reported here), the Doppler width is smaller than the 
Stark width and Doppler convolution is no longer necessary. 

4. Conclusions 

This study shows that Stark broadened infrared Brackett lines of hydrogen follow the same 
trends as the lines of lower series which are more known theoretically and experimentally. 
This study will allow missing absorption in the spectra of stellar atmospheres in the infrared 
due to the lack of data to be filled in a next future. Dedicated experiments and comparisons 
with other methods, for instance FFM [32], would be helpful to test these theoretical results, 
which will be also constrained by observational data. The corresponding tables will be 
accessible at |http://amrel.obspm.fr/stark-h 
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Figure 4: MMM Br(3 line at 10 15 cm 3 , versus Act = AX/F without Doppler convolution 
for different temperatures (red : 1000 K; blue with empty circles : 2500 K; violet with empty 
squares: 10000K; green: 19950K). 



9 




Figure 5: MMM Bra line at N e =10 12 cm 3 ; with (red) and without (blue) Doppler broadening, 
and the asymptotic Holstmark limit (black) in |Aa|~ 5//2 for 11 different temperatures in K 
(1000, 2500, 5000, Iff, 2 10 i , 3.98 10*, 7.94 1-58 10 5 , 3.16 10 5 , 6.31 10 5 , 1.26 l(fi). 
The corresponding values of the profiles in the line center are respectively : 0.28 ; 0.30 ; 0.34 
; 0.46 ; 0.56 ; 0.69 ; 0.88 ; 1.23 ; 1.45 ; 1.87, without Doppler effect (blue), and : 0.102 ; 
0.071; 0.053; 0.039; 0.028 ; 0.020; 0.014; 0.010 ; 0.0072 ; 0.0051; 0.0036 including Doppler 
effect (red). 
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Figure 6: Same as Fig\^ but for Brf3, The values of the profiles in the line center are, for 
the 11 different temperatures of Fig\5\ respectively : 0.124 ; 0.129 ; 0.136 ; 0.149 ; 0.168 ; 
0.195 ; 0.233; 0.285 ; 0.355 ; 0.451 ; 0.579 without Doppler effect (blue), and 0.099; 0.082 
; 0.067 ; 0.053 ; 0.040 ; 0.030 ; 0.0215 ; 0.015 ; 0.011 ; 0.0078; 0.0056 including Doppler 
effect (red). 
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Figure 7: Same as Fig!^5\ but for Bry, The values of the profiles in the line center are, f for 
the 11 different temperatures of Fig^ respectively : 0.107 ; 0.107 ; 0.109 ; 0.113 ; 0.121 ; 
0.132 ; 0.150 ; 0.175 ; 0.211 ; 0.260 ; 0.327 without Doppler effect (blue), and 0.091 ; 0.080 
; 0.069 ; 0.057 ; 0.0450 ; 0.034 / 0.025 ; 0.018 ; 0.013 ; 0.0094 ; 0.0067 including Doppler 
effect (red). 
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